arXiv:math/0211344vl [math.OC] 21 Nov 2002 


ALGEBRAIC METHODS FOR COMPUTING SMALLEST ENCLOSING 
AND CIRCUMSCRIBING CYLINDERS OF SIMPLICES 

RENE BRANDENBERG AND THORSTEN THEOBALD 


Abstract. We provide an algebraic framework to compute smallest enclosing and small¬ 
est circumscribing cylinders of simplices in Euclidean space E". Explicitly, the compu¬ 
tation of a smallest enclosing cylinder in E^ is reduced to the computation of a smallest 
circumscribing cylinder. We improve existing polynomial formulations to compute the 
locally extreme circumscribing cylinders in E^ and exhibit subclasses of simplices where 
the algebraic degrees can be further reduced. Moreover, we generalize these efficient for¬ 
mulations to the n-dimensional case and provide bounds on the number of local extrema. 
Using elementary invariant theory, we prove structural results on the direction vectors 
of any locally extreme circumscribing cylinder for regular simplices. 


1. Introduction 


Radii (of various types) belong to the most important functionals of polytopes and 
general convex bodies in Euclidean space E” [S, |T5,0 , and they are related to applications 


in computer vision, robotics, computational biology, functional analysis, and statistics 
T5[|). Following the notation in Q, the outer j-radius Rj{V) of a convex body 


see 


C C E” is the radius of a smallest enclosing j-dimensional sphere in the optimal orthogonal 
projection of C onto a j-dimensional linear subspace. Studying these radii, mainly for 
regular simplices and regular polytopes, is a classical topic of convex geometry (see |^, 

00 ). 


From the computational point of view, most of the existing algorithms for computing 
these radii focus on approximation A major reason is that exact computations 

lead to algebraic problems of high degree, even for computing, say, the outer (n—l)-radius 
in E" (already if n = 3). However, since some approaches for computing radii of general 
polytopes consider the computation of a smallest enclosing cylinder of a simplex as a 
black box within a larger computation B. 0. these core problems on simplices are of 
fundamental importance. 

Recently, the authors of demonstrated that using their state-of-the-art numerical 
polynomial solvers, various problems related to cylinders in E^ can be solved rather effi¬ 
ciently. In particular, the authors give a polynomial formulation for the smallest circum¬ 
scribing cylinder of a simplex in E^, whose Bezout number - the product of the degrees 
of the polynomial equations - is 60. However, these equations contain certain undesired 
solutions with multiplicity 4, and as a result of these multiplicities the computation times 
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(using state-of-the-art numerical techniques) are about a factor 100 larger than those of 
similar problems in which all solutions occur with multiplicity 1. 

Here, we provide a general algebraic framework for computing smallest enclosing and 
circumscribing cylinders of simplices in E”. First we reduce the computation of a smallest 
enclosing cylinder in to the computation of a smallest circumscribing cylinder, thus 
combining these two problems. Then we investigate smallest circumscribing cylinders of 
simplices in E^. We improve the results of [Q by providing a polynomial formulation for 
the locally extreme cylinders, whose Bezout bound is 36 and whose solutions generically 
have multiplicity one. Our formulations use techniques from the paper which studies 
the lines simultaneously tangent to four unit spheres. These techniques also enable us 
to present classes of simplices for which the algebraic degrees in computing the smallest 
circumscribing cylinder can be considerably reduced. 

Then, in Section H, we give a generalization of our approach to smallest circumscribing 
cylinders of a simplex in E”. Based on this formulation we give bounds on the number 
of locally extreme cylinders based on the Bezout number. Since this bound is not tight, 
we provide better bounds for small dimensions; these bounds are based on mixed volume 
computations and Bernstein’s Theorem. Moreover, we study in detail the locally extreme 
circumscribing cylinders of a regular simplex in E”. To exploit many symmetries in the 
analysis, we provide a formulation based on symmetric polynomials. Using elementary 
invariant theory we show that the direction vector of every locally extreme circumscribing 
cylinder has at most three distinct values in its components. With this result we can 
illustrate our combinatorial results on the number of solutions for general simplices. 

As a byproduct of our computational studies, we discovered a subtle but severe mistake 
in the paper on the explicit determination of the outer (n—l)-radius for a regular 
simplex in E"^, thus completely invalidating the proof given there. In the appendix we 
give a description of that flaw, including some computer-algebraic calculations illustrating 
it. 


2. Preliminaries and background 


2.1. j-radii and cylinders. Throughout the paper we work in Euclidean space E”, i.e., 
M" with the usual scalar product x ■ y = norm ||x|| = {x ■ We write 

x"^ for X ■ X. 

A j-flat is an affine subspace of dimension j. For a convex polytope P c E"" (or a hnite 
point set P C E"") and a j-flat E, we consider 

TZV{V,E) := maxdist(p, E), 

p£P 


where dist(p, E) denotes the Euclidean distance from p to E. The outer j-radius of V is 


nv{v,E). 

tj IS an (n—j)-nat 

The choice of the indexing in the j-radius stems from the fact that it measures the 
radius of an enclosing j-dimensional sphere in the optimal orthogonal projection of V 
onto a j-dimensional linear subspace (cf. iHl)- 
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One of the most natural representatives of this class is the one with j = 2, n = 3, i.e., 
the smallest enclosing (circular) cylinder of a polytope. In E", we define a cylinder to be 
a set of the form 

bd(£ + pB^), 

where £ is a line in E**, denotes the unit ball, p > 0, the addition denotes the Minkowski 
sum, and bd(-) denotes the boundary of a set. We say that P can be enclosed in a cylinder 
C if P is contained in the convex hull of C. Thus the outer (n—l)-radius gives the radius 
of the smallest enclosing cylinder of a polytope. 

A simplex in E”^ is the convex hull of n + 1 affinely independent points. An enclosing 
cylinder C of a simplex V is called a circumscribing cylinder of V if all the vertices of V 
are contained in (the hypersurface) C. 

2.2. Smallest circumscribing cylinders and smallest enclosing cylinders. The fol¬ 
lowing statement connects the computation of a smallest enclosing cylinder of a polytope 
with the computation of a smallest circumscribing cylinder of a simplex.^ 

Theorem 1. Let V = {pi,... ,Pm} be a set of m > 4 points in E^, not all collinear. IfV 
can be enclosed in a circular cylinder C of radius r, then there exists a circular cylinder 
C of radius r enclosing all elements ofV such that the surface C passes through 

(i) at least four non-collinear points ofV, or 

(ii) three non-collinear points ofV, and the axis i of C is contained in 

(a) the cylinder naturally defined by spheres of radius r centered at two of these 
points; 

(b) the double cone naturally defined by spheres of radius r centered at two of 
these points (and these spheres are disjoint); 

(c) or the set of lines which are tangent to the two spheres of radius r centered at 
two these points and which are contained in the plane eguidistant from these 
points (and the spheres are non-disjoint). 

Moreover, C can be transformed into C by a continuous motion. 

Figures |l| and || visualize the three geometric properties in the second possibility. 

Since the second possibility in Theorem |l| characterizes the possible special cases, this 
lemma in particular reduces the computation of a smallest enclosing cylinder of a simplex 
in E^ to the computation of a smallest circumscribing cylinder of a simplex. Namely, 
it suffices to compute the smallest circumscribing cylinder (corresponding to case (i)) as 
well as the smallest enclosing cylinders whose axes satishes one of the condition in (ii); 
the latter case gives a constant number of problems of smaller algebraic degree (since the 
positions of the axes are very restricted). 

Remark 2. Before we start with the proof, we remark that Theorem ^ and its different 
cases show a quite similar behaviour as the well known statement that the (unique) 
circumsphere of a simplex touches all its vertices, or one of its great (n—l)-circles is the 
circumsphere of one of the (n—l)-faces of the simplex (see p. 54]). 

^We remark that a similar statement has already been used in |^, but the manuscript referenced 
there does not contain a complete proof. 
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Figure 1. Extreme situations of the set of hyperboloids for disjoint spheres 




(b) Degenerated hyperboloid for Xh = 


(a) Hyperboloid for 0 < < 2r^/a 

Figure 2. The left hgure shows a general situation for disjoint spheres; 
the right hgure shows an extreme situation for non-disjoint spheres 


a/2 


In the proof we will apply the following geometric equivalence. A point a; G is 
enclosed in a cylinder with axis £ if and only if £ is a transversal of the sphere with radius 
r centered at x (i.e., £ is a line intersecting the sphere). 

Proof of Theorem |^. Let C be a cylinder with axis i and radius r enclosing V. Then, 
denoting by Si := S{pi, r) the sphere with radius r centered at Pi, £ is a common transversal 
to Si,Sm- By continuously translating and rotating i, we can assume that i is tangent 
to two of the spheres, say Si and 5*2. Further, by changing coordinates, we can assume 
that Si and S 2 have the form Si = *S'((0, 0, 0)^, r), S 2 = S{{a, 0, 0)^, r) for some a > 0. 

The set of lines tangent to two spheres of radius r constitutes a set of hyperboloids 
(see, e.g., iH)- Moreover, any of these hyperboloids touches the sphere Si on a circle 
lying in a hyperplane parallel to the i/z-plane. Hence, the set of hyperboloids can be 
parametrized by the x-coordinate of this hyperplane which we denote by Xh- 

If S'! n 5*2 = 0 then the boundary values are Xh = 0 and Xh = ‘lr‘^ ja. These two 
extreme situations yield a cylinder and a double cone with apex (a/2, 0,0)^, respectively 
(see Figure |l[). For 0 < Xh < ‘lr‘^ ja we obtain a hyperboloid of one sheet (see Figure ia)). 



























SMALLEST ENCLOSING AND CIRCUMSCRIBING CYLINDERS 


5 


If S'inS '2 7 ^ 0 then the boundary values are Xh = ^ and Xh = a/2. Here, for 0 < Xh < a/2 
we obtain hyperboloids of one sheet, too. For Xh = a/2 the hyperboloid degenerates to 
a set of tangents which are tangents to the circle with radius = a/ 4r^ — a? in the 
hyperplane x = a/2 (see Figure |Kb)). 

Let Xh^Q be the parameter value of the hyperboloid containing the line The tangent to 
Si and S2 is contained in the hyperboloid with some parameter value Xhfi- By decreasing 
the parameter Xh starting from Xh^ the hyperboloid changes its shape towards the cylinder 
around Si and 5 * 2 . Let Xh^i be the infimum of all 0 < < Xh,o such that the hyperboloid 

does not contain a generating line tangent to some other sphere S{pi, r) for some 3 < i < 
m. If Xh,i = 0, then by choosing any point of V not collinear to pi and p 2 we are in case 
(ii) (a). 

If Xh,i > 0 then let ps be the corresponding point. Let T(S'i, 52 , S' 3 ) denote the set 
of lines simultaneously tangent to 5i, S 2 , and S^. Now let Xh .,2 be the infimum of all 
0 < T/i < Xhp such that there exists a continuous function £ : {xh^ 2 -i Xh,i) T{{Si, S 2 , Ss}) 
with £{xh) lying on the hyperboloid with parameter Xh- Since the spheres are compact, 
the inhmum is a minimum. If Xh ,2 > 0 then one of three hyperboloids involved by the 
three pairs of spheres must be one of the extreme hyperboloids in that situation and we 
are in cases (ii) (a), (b), or (c). If Xh ,2 = 0 then we distinguish between two possibilities. 
Either during this process we also reached a tangent to some other sphere S{pi,r) for 
some 4 < z < m; in this case we are in case (i). Or during the transformation all the 
points P 4 ,... ,Pm are enclosed in the cylinder with axis i and radius r, but none of them 
is contained in it. Then we arrive at situation (ii) (a). □ 

3. Computing the smallest circumscribing cylinders of a simplex in 

So far, we have seen how to reduce the computation of a smallest enclosing cylinder of a 
simplex in E^ to the computation of a smallest circumscribing cylinder. In order to apply 
algebraic methods to compute a smallest circumscribing cylinder, there are many different 
ways to formulate that problem in terms of polynomial equations. It is well-known that 
the computational costs of solving a system of polynomial equations are mainly dominated 
by the Bezout number (= product of the degrees) and the mixed volume (the latter one 
is discussed in Section p. See BBS for comprehensive introductions and the state- 
of-the-art. Hence, it is an essential task to find the right formulations. Moreover, we are 
interested in simplex classes for which the degrees can be further reduced. 

3.1. General simplices in E^. In the proof of P, Theorem 6], a polynomial formulation 
is given to compute the smallest enclosing cylinder of a simplex in E^. This formulation 
describes the problem by three equations in the direction vector v = {vi,V 2 ,V 3 )'^ of the 
line, one of them normalizing the direction vector v by 

(3.1) vl + vl + vl = 1. 

The equations are of degree 10, 3, and 2, respectively, thus giving a Bezout number of 
60. However, as pointed out in that paper, some of the solutions to that system are 
artificially introduced by the formulation and occur with higher multiplicity, and there 
are only 18 really different solutions. Even more severely, in the experiments in that paper 



6 


R. BRANDENBERG AND T. THEOBALD 


(using Synaps, a state-of-the-art software for numerical polynomial computations), the 
numerical treatment of these multiple solutions needs much time, roughly a factor 100 
compared to similar systems without multiple solutions. 

Here, we present an approach, which reflects the true algebraic bound of 18. Namely, we 
give a polynomial formulation with Bezout bound 36 in which every solution generically 
has multiplicity one. The additional factor 2 just results from the fact that due to the 
normalization condition (|3.1|) every solution v also implies that —v is a solution as well. 

Our framework is based on in which the lines simultaneously tangent to four unit 
spheres are studied. A line in is represented by a point u G E^ lying on the line and a 
direction vector u G E^ with = 1. We can make u unique by requiring that u • u = 0. 
A line £ = (u, v) has Euclidean distance r from a point p G E^ if and only if the quadratic 
equation [u + tv — pY = has a solution of multiplicity two. This gives the condition 


{v ■ {u-p)y 


— {u — pY +r‘^ = 0 . 


Expanding this equation yields 


(3,2) 


2„,2 


V U 


— 2v u ■ p + V p — {v ■ p) — r V = 0 . 


Rather than using = 1 to further simplify this equation, we prefer to keep the homoge¬ 
nous form, in which all terms are of degree 4. 

Now let pi,... ,p 4 be the affinely independent vertices of the given simplex. Without 
loss of generality we can choose p 4 to be the located in the origin. Then the remaining 
points span E^. Subtracting the equation for the point in the origin from the equations for 
Pi) P 2 ) P 3 gives the following program to compute the square of the radius of the minimal 
circumscribing cylinder. 


(3.3) 


min 

S.t. U ■ V 

2v‘^u ■ Pi 


0 , 

v‘^p‘f — {v ■ piY , 1 < i < 3 , 

1 . 


We remark that the set of admissible solutions is nonempty; a proof of that 
(for general dimension) is contained in Section 

Since the points pi,P 2 ,P 3 are linearly independent, the matrix M := (pi, 
invertible, and we can solve the equations in the bottom line of ( p.3|) for u: 


statement 
P‘2,Pzf is 


(3.4) 


u 



v‘^p‘1 — {v ■ PiY \ 

uVi - ■ P2Y 

V^pI - {v ■ P3)2 J 


Now substitute this expression for u into the the objective function and into the hrst 
constraint of the system (|3.3|) . After setting = 1 in the denominator of the hrst 
constraint, this gives a homogeneous cubic equation which we denote by gi{vi,V 2 , U 3 ) = 0 . 
Hence, we arrive at the following polynomial optimization formulation in terms of the 
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variables V 2 i and v^. 

/ / \Y 

min ^ v‘^pI — {v ■ P 2 Y 

(3.5) y y v'^pI - (n • pY j j 

s.t. gi{vi,V 2 ,v:i) = 0, 

g2{vi,V2,v^) ■.= -I = 0 . 

Note that the objective fnnction is a homogeneons polynomial of degree 4. We denote 
this polynomial by /. 

Using Lagrange mnltipliers Ai and A 2 , a necessary local optimality condition is 


(3.6) grad / = Aigrad gi + Aagrad g 2 . 

By thinking of an additional factor Aq before grad / and considering (|3.6|) as a system of 
linear eqnations in Aq, Ai, A 2 , we see that if (|3.6|) is satisfied for some vector v then the 
determinant 

df dgi dg2 \ 
dvi dvi dvi ' 

df dgi dg2 
dv2 dv2 dv2 
df dgi dg2 
dv3 dv3 dv3 

vanishes. 

Lemma 3. (a) Any direction vector {v\,V 2 ,v^Y ^ o,xis of a locally extreme 

circumscribing cylinder satisfies the polynomial system (p3|) . If this system has only 
finitely many solutions then this number is bounded by 36. 

(b) For a generic simplex the number of solutions is indeed finite, and all solutions have 
multiplicity one. 


(3.7) 


- 


det 


Proof. Let v be the direction vector of an axis of a locally extreme circnmscribing cylinder. 
Then v satisfies the hrst constraint of and the determinant (E3) vanishes. Since 

these are homogeneons eqnations of degree 3 and 6, respectively, Bezont’s Theorem implies 
that in connection with = 1 we obtain at most 36 isolated solntions. 

For the second statement it snffices to check that for one specihc simplex there are only 
hnitely many solntions and that all solntions are pairwise distinct. □ 


3.2. Special simplex classes in E^. In this section, we investigate conditions nnder 
which the degree of the resnlting eqnations is rednced. Moreover, we show that for the 
eqnifacial simplex, the minimal circnmscribing radins can be compnted qnite easily. 


We nse the following classihcation from |21, 22 


Proposition 4. Let T be a simplex in E^ with vertices pi,... ,p 4 . The polynomialgi in the 
cubic equation factors into a linear polynomial and an irreducible quadratic polynomial if 
and only if the four faces ofT can be partitioned into two pairs of faces {Fi, F 2 }, {F 3 , F 4 } 
with area{Fi) = area{F 2 ) 7 ^ area{Ffi) = area^Ffi). Moreover, gi factors into three linear 
terms if and only if the areas of all four faces ofT are equal. 
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First let us consider the case where gi decomposes into a linear polynomial and an irre¬ 
ducible quadratic polynomial. By optimizing separately over the linear and the quadratic 
constraint, the degrees of our equations are smaller than for the general case. Namely, 
analogously to the derivation in Section |3.1j , for the quadratic constraint we obtain a 
Bezout bound of 

(3 + 1 + 1) ■ 2 ■ 2 = 20 , 
and for the linear constraint we obtain 


(3 + 0 + 1)-1-2 = 8, 


Thus, we can conclude: 


Lemma 5. If the four faces of the simplex can he partitioned into two pairs of faces 
{Fi,F 2 }, {+ 3 ,+ 4 } with area{Fi) = area{F 2 ) 7 ^ area{F^) = area^F^) then there are at most 
28 isolated local extrema for the mimimal circumscribing cylinder. They can be computed 
from two polynomial systems with Bezout numbers 20 and 8 , respectively. 


Equifacial simplices. A simplex in is called equifacial if all four faces have the same 
area. By Proposition for an equifacial simplex the cubic polynomial gi factors into three 
linear terms. Hence, we obtain at most 3 ■ 8 = 24 local extrema. Somewhat surprisingly, 
using a characterization from it is even possible to compute smallest circumscribing 
cylinder of an equifacial simplex esentially without any algebraic computation. 

Namely, it is well-known that the vertices of an equifacial simplex T can be regarded 
as four pairwise non-adjacent vertices of a rectangular box (see, e.g., fl^). Hence, there 
exists a representation pi = {wi,W 2 ,W 3 )'^, p 2 = {wi, 

P 4 = {-Wi, -W 2 , Ws)'^ with Wi, W 2 , W 3 > 0. 


-W2,-W3)^, P3 = {-Wi,W2,-W3f 


Assuming without loss of generality = 1, 


gives 


(3.8) 


{v ■ Pif + 2 u- Pi = ^2 w] + 1 < * < 4 . 

Subtracting these equations pairwise gives 

A{w2U2 + W3U3) = -A{wiW3ViV3 + W1W2V1V2) 

(for indices 1, 2) and analogous equations, so that 

WiUi = -W2W3V2V3, W2U2 = -W1W3V1V3, W3U3 = -W1W2V1V2. 

Since n ■ n = 0, this yields V\V 2 V 3 = 0. Without loss of generality we can assume vi = 0. 
In this case, 

W 2 W 3 


u = 


Wi 


-^ 2 ^ 3 , 0,0 


So we can express (|3.8|) in terms of the direction vector n. 


2 2 I 22 2 I 

W2V2 + = 2 ^Wj + 


W 2 W 3 \ 2 

- V 2 V 3 I - r 
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which, by using V 2 + u 


'I 


1 , gives 


(3.9) 


r 


,2 



2 I 2 I 2 
V 2 + + W2- 


Thus, by computing the derivative of this expression = r^(u 2 ) and taking into account 
the three cases u, = 0, we can reduce the computation of the minimal circumscribing 
cylinder to solving three univariate equations of degree 3. However, we can still do better. 
Substitute Z 2 ■= V 2 , and let p be the expression for in terms of 2 ) 2 , 



Since the second derivative of that quadratic function is negative, ^(^ 2 ) is a concave 
function. Hence, within the interval Z 2 € [0,1], the minimum is attained at one of the 
boundary values Z 2 G {0,1}. Consequently, two of the components of (ui,U 2 ,U 3 )^ must 
be zero and therefore v is perpendicular to two opposite edges. Since the latter geometric 
characterization is independent of our specific choice of coordinates, we can conclude: 

Lemma 6. If all four faces of the simplex T have the same area then the axis of a 
minimum circumscribing cylinder is perpendicular to two opposite edges. 

Hence, for an equifacial simplex it suffices to investigate the cross products of the three 
pairs of opposite edges (equipped with an orientation), and we do not need to solve a 
system of polynomial equations at all. 

In order to illustrate how these three solutions relate to the 18 solutions of the general 
approach above, we consider the regular simplex in E^. In the general approach, as already 
pointed out in |^, the six edge directions piPj (1 < i < j < 4) all have multiplicity 1, 
and each of the three directions in Lemma |^, pip2 x P 3 P 4 , pip^ x P2P a, PiPa x P2P3, have 
multiplicity 4. 

4. Smallest circumscribing cylinders in higher dimensions 

In Section we have given polynomial formulations with small Bezout number for com¬ 
puting smallest circumscribing cylinders of a simplex in E^. Using the characterization 
in of lines simultaneously tangent to 2n—2 spheres in E", we generalize these formu¬ 
lations to smallest circumscribing cylinders of a simplex in E”, n > 2 . Analogous to the 
three-dimensional case let pi, ■ ■ ■ ,Pn+i be the affinely independent vertices of the simplex 
in E”, and let Pn+i be located in the origin. 

First note that (|3.3|) also holds in general dimension n if we replace the index 3 by the in¬ 
dex n. Since the points pi,... are linearly independent, the matrix M := (pi,... ,Pn)'^ 
is invertible, and we can solve for u\ 


(4.1) 
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Hence, by generalizing the formulation for the three-dimensional case, we obtain the 
program 


— {v ■ 

vVn-iv-Pn? 
s.t. gi{vi,...,Vn) = 0, 

g2{vi, ... ,Vn) := - 1 = 0, 

where gi denotes the cubic equation as before. In order to show that set of admissible 
solutions for our optimization problem is nonempty, we record the following result. 

Lemma 7. For any simplex in E"’ the (”2 edge directions of the simplex are direction 
vectors of circumscribing cylinders. 

Proof. Since the edge directions pi — pj have a simple description in the basis pi, ■ ■ ■ ,Pn, 
we express the cubic equation gi{v) = 0 in that basis. Let v be an arbitrary direction 
vector, and let the representation of v in the basis pi,... ,pn be 

n 

V = tiPi. 

i=l 




Further, let p[,...,p'^ be a dual basis to pi,...,pn] i.e., let p[,...,p'.^ be dehned by 
p'i ■ pj = 6ij, where 6ij denotes Kronecker’s delta function. By elementary linear algebra, 
we have U = p[- v. 

When expressing u in this dual basis, u = J2KPii second constraint of (|3.3|) gives 

' 1/22 / \ 2 \ 

“i = ^ ) • 

Substituting this representation of u into the equation gi{v) = 0 gives 


0 = gi{v) = v^{u-v) = 


■ V = v 


5 ^ u'iU , 


, 2=1 


2=1 


where the last step uses the duality of the bases. Hence, we obtain the cubic equation 

1 


- {v ■ pif)ti = 0. 


2=1 


Expressing v in terms of the t-variables yields 

1 


^ ^ ^ ^ fdijktitjt}^ 0 , 




l<i<j <k<n 
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where 


CXij 


Pijk 


{vo\2{Pi,Pj)Y 


f Pi -Pi Pi- Pj 

V Pj ■ Pi Pj ■ Pj 




det 


Pi ■ Pj Pi ■ Pk 
Pk ■ Pj Pk ■ Pk 


+ det 


Pi -Pk Pi- Pj \ 
Pj - Pk Pj - Pj J 


+ det 


Pj - Pk Pj - Pi 
Pi -Pk Pi- Pi 




and vo\ 2 {pi,Pj) denotes the oriented area of the parallelogram spanned by Pi and pj. In 
terms of the t-coordinates, the (” 2 ^) edges of the simplex are t = e*, 1 < i < n, and 
t = Ci — Cj, I < i < j < n, where e* denotes the i-th standard nnit vector. For all these 
edges, the cnbic eqnation is satisfied. □ 


Considering Lagrange mnltipliers Ai and A 2 yields the following necessary optimality 
condition. 


(4.3) 


grad / 

g2(vi, ...,Vn) 


Aigrad gi + A 2 grad g 2 , 

0 , 

0 . 


Since the Bezont bonnd of this system is 3” ■ 3 ■ 2 = 2 ■ we have: 


Lemma 8. Forn > 2, the number of isolated local extrema for the minimal circumscribing 
cylinder is bounded by 2 ■ . 


This bonnd is not tight. Trying to rednce this npper bonnd of isolated solntions like 
in the three-dimensional case, we can eliminate the linear occurrences of the Lagrange 
variables Ai and A 2 . Generalizing (|3.7|) , we have to consider the vanishing of all 3x3- 
subdeterminants of the matrix 


(4.4) 


/ 

dvi 

-M. 

dv2 


\ dVn 


dgi 

dg 2 

dvi 

dvi 

dgi 

dg 2 

dv 2 

dv 2 

dgi 

dg 2 

dVn 

dVn 


Thus, for n > 4 we arrive at a non-complete intersection of equations where we have more 
equations than variables. Hence, we cannot apply our Bezont bound on these systems. 

However, for small dimensions we can improve Lemma ^ by directly working on the 
formulation (|4.3|) . In order to provide better bounds, we use well-known characterizations 
of the number of zeroes of a polynomial equation by the mixed volume of a Minkowski 
sum of polytopes (for an easily accessible introduction into this topic we refer to 0 ]). 
Here, let C* := C\ {0}. 


Lemma 9. For 2 <n <7, the number of solutions of the system (O) in (ui,..., Ai, A 2 ) 
e (C*)”+^ is bounded by 


n -I- 1 
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where {”} denotes the Stirling number of the second kind (see, e.g., [^, [2^ 
The sequence starts as follows 


n 

2 

3 

4 

5 

6 

7 

NTT 

6 

36 

150 

540 

1806 

5796 


Proof. For a polynomial h = G C[xi,..., x^], let 

NP(/i) := conv{a G Mg : 7 ^ 0} 


denote the Newton polytope of h (see, e.g., 0, §7.1]). Let hi,... ,hn be the polynomials 
of the gradient equation in ([4.3|) . Further let Pi,..., P„, Qi, Q 2 be the Newton polytopes 
of hi,..., hn, gi, g 2 for generic instances of these equations. 

Recall that the mixed volume MV(Pi,..., P„, Qi, Q 2 ) is the coefficient of the monomial 
ArA 2 • ■ ■ An-/ir/i 2 in the (n+ 2 )-dimensional volume Vol„+ 2 (AiPi + .. ■+XnPn+PiQi+P 2 Q 2 ) 
(which is a polynomial expression in Ai,..., A„, pi,/i 2 ). By Bernstein’s Theorem, the 
number of isolated common zeroes in (C*)"’*'^ of the set of polynomials hi,..., hn,gi,g 2 
is bounded aboved by 

MV(Pi,...,P„,Qi,Q2) 


(see 10, Chapter 8 , Theorem 5.4]). For every given n this volume can be computed using 
software for computing mixed volumes (see, e.g. 


12 , 22 


□ 


We conjecture that for any n> 2, the number of isolated solutions in (C*)"'”''^ is bounded 
by 6 {”y}. 


4.1. The regular simplex in E". Here, we analyze the local extrema of circumscribing 
cylinders for the regular simplex. Our aim is both to illustrate the algebraic formulations 
given before and to relate our investigations to classical investigations on the regular 
simplex in convex geometry. In order to achieve many symmetries in the algebraic for¬ 
mulation, we use a slightly modified coordinate system that is particularly suited for the 
regular simplex; these coordinates have also been used in |^, 0 ] . 

The equation xi -|- ... -f- Xn+i = 1 defines an n-dimensional affine subspace in 
Now let the regular simplex in this n-dimensional subspace be given by the n -|- 1 vertices 
Pi = Ci, where e* denotes the nth standard unit vector, 1 < z < n -|- 1. We consider the 
tangency equation (|3.2|) for the point Pn+i, 

22o2 i2 2 22 

V u — 2v Un+i + V — — r V = U . 

Subtracting this equation from the equation for pi, 1 < z < n, yields 

2v'^{ui - Un+i) = -{v} - u^+i), 1 < z < n . 

Moreover, the embedding into the hyperplane x, = 1 implies Ui = 1. In order 

to solve these n -f 1 equations for u, let M be the (n -f 1) x (n -f l)-matrix whose z-th row 
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contains the vector ef — and whose n-th row is (1,1,..., 1). Since M is invertible, 
we obtain 


(4.5) 




/ -(yl - v^+i) \ 

V 


As before, substituting this expression into u-v = 0 and setting = 1 in the denominator 
gives a cubic equation gi{v) = 0. Hence, we obtain the following optimization problem. 
Here, the objective function / stems from the condition for the vertex Pn+i, and the 
condition n* = 0 comes from the embedding. 


(4.6) 


s.t. 


min - 2un+i + 1 - 

gi{vi,...,Vn+i) = 0 , 

n+l 

= 0 > 

2=1 

= 1 . 


First we record that the functions / and gi are symmetric polynomials in the variables 
Ti,..., Vn+i- In order to show this, let ai,..., (Jn+i be the elementary symmetric functions 
in vi,.. .,Vn+i, 


(Ti = Ti + . . . + Vn+l , 

^ ■ ■ ■ Tjj, , 

l<2l<...<2fc<r2H-l 

^n+l ' ‘ * ”^22+1 


(see, e.g., § 0 ). By providing explicit expressions for / and gi as polynomials in the 
elementary symmetric polynomials ui,..., cr„+i, the symmetry of / and gi follows. More 
precisely, we obtain: 


Lemma 10. The quartic polynomial /(ui,..., Un+i) and the cubic polynomial gi{vi,..., 
n„+i) are symmetric polynomials in the variables Vi,, Vn+i- In terms of the elementary 
symmetric functions, f results in 


f 


- ^nala 2 + 2{n - l)al 


Aal + 8(72 + 4n) + cxicra - 0-4 , 


and the homogeneous polynomial gi results in 

1 


9i = 


2 (n+ 1 ) 


(-(n - 2)af + 3{n - l)aia 2 ) - - 0-3 . 
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Since di = 0 and ~ 2 ct 2 , we can also deduce the following formulation of 

our optimization problem: 

Corollary 11. Finding the critical values of the minimization problem ( |4.6|) is equivalent 
to finding the critical values (ui,..., Vn+i)'^ of the maximization problem 


max (74 

S.t. (Ti = 0, 

(4-7) as = 

as = 0, 

where ai are the elementary symmetric functions in vi,, Un+i- 

Theorem 12. The direction vector (ui,... ,Vn+i)'^ of any locally extreme circumscribing 
cylinder satisfies |{ni,... ,r;n+i}| < 3, i.e., for each solution vector the components take 
at most three distinct values. 


Proof. For n < 2, the statement is trivial, so we can assume n > 3. Let v be the direction 
vector of a locally extreme circumscribing cylinder with = 1. Using Corollary 0 let 
f(v) := -afiv), gi{v) := afiv), g 2 {v) := a 2 {v) - 1 / 2 , and 5 ( 3 ( 1 ;) := ai(i;). As a necessary 
condition for a local extremum, for any pairwise different indices a,b,c,d G { 1 ,..., n + 1 } 
the determinant 



/ 

dVa 

dgi 

dg2 

dgs \ 

dVa 


dVa 

dVa 



dgi 

dg2 

dga 

(4.8) det 

dvi, 

dvf, 

dgi 

dvh 

dg2 

dvh 

dg3 


dvc 

dvc 

dvc 

dvc 


[ -PL 

\ dva 

dgi 

dg2 

dgs , 
dvd ' 


dvd 

dvd 

vanishes. Since /, gi, 32 , and 5(3 

are symmetric functions 

we can assume without loss of generality 

a = 

1 , 

6 = 2 , 

«n:=Er=^iand/3. = ESu^ 

we can write 




and d = 4. Setting 


dgs 

dvi 

dg 2 

dvi 


dgi 

dvi 


1 , 


4 

jV* 

^ VjVk + an^Vj + - [al - (3n) 

l<j<k<A 3=1 


(1 < z < 4). Moreover, since a 3 (n) = 0, we can consider as + ^ instead of This 
allows to express the resulting expression easily in terms of and fin- More precisely, we 
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obtain 


df 

a 3 + — = Vi 
OVi 


/ 


\ 


VjVk + anY^j+o 


1 l<j<fc<4 


J = 1 




Thus we can consider the determinant (|4.8|) as a polynomial in fi, ^ 2 , T 3 , ^ 4 , Ofn, /^n- Evaluat¬ 
ing this 4 X 4-determinant A shows that it is independent of an, /5„ and that it factors 
as 

A = (m - V 2 ){vi - V3){vi - V^){V 2 - V3){V2 - Vi){v3 - V^) . 

Hence, |{ni, V 2 , V 3 , n 4 }| < 3, and this holds true for any quadruple (a, b, c, d) of indices. □ 


Using this result, we illustrate the occurrence of the Stirling numbers in Lemma ^ for the 
case of a regular simplex. There are {"' 3 ^} ways to partition the set U := {ui,... 
into three nonempty subsets Ui, V 2 , V 3 . We assume that Vi E Vi, 1 < i < 3, and that all 
variables within the same set take the same value. Setting k := |l/i| and I := IV 2 I, the 
formulation in Corollary yields the system of equations 


(4.9) 


kvi + lv 2 + {n + 1 - 
kvf + lvl + {n + 1 — 


0<ii<i2<'i3<3 

n+»2+®3=3 


H 


l \ kn + 1 — k — I 
*2/ V U 


k - 1 )V 3 

k — l)vl 


0 , 

1 , 

0 . 


If one of the indices fc, /, or n -|- 1 — fc — / is zero then this system consists of three 
equations in two variables, so we do not expect any solutions. For every choice of k, I 
corresponding to a partition into nonempty subsets, we obtain a system of equations with 
Bezout number 6 . Thus, whenever the values of vi, V 2 , and V 3 in the solutions to ( [4.9|) 
are distinct, then this reflects the bound in Lemma 

In particular, in the case n = 4 we obtain the following 150 solutions. 


k = 1, I = 1: The six solutions for (vi, ^ 2 , ^ 3 ) of the system (|4.9|) are 

T 


1 1 


— V 110 -30*^15,—V 110 +30i^,- — 715) , 
20 ^ ’ 20 ^ ’ 10 / ’ 


and the solutions obtained by permuting the hrst two components of the hrst 
solution and by changing the signs and/or permuting the hrst two components in 
the second solution. 

For the program ([4.7|) in the variables (ui,..., us)^, this gives ( 2 ) (^) = 20 critical 
positions of the form (i.e., up to variable permutations) 

T 

:, 0 , 0 , 0 ' 


1 

72 ’ 


1 

72 ’ 


20 complex solutions of the form 


(^-^^110-30*715 ,-^^110 +30*715, ^715, ^715, 
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and 20 complex solutions of the form 
( —V 110 - 30iv^, —V 110 + 30iv^, v^, 

V 20 V ’20^ ’ 10 ’ 10 ’ 10 

A; = 1, / = 2: Here, we obtain 30 solutions of the form 



/II 1 i\^ 

O’2' 2 ’ 2' 2 J • 


30 solutions of the form 









and 30 solutions of the form 






The global mininum is attained for the vector (O, and the objective value 

of the global optimum is 49/80. Hence, the radius of the smallest circumscribing cylinder 
for a regular simplex in with edge length a/2 is ^49/80 = 7v/5/20 ^ 0.7826 . 


Appendix: An error in the results of Weissbach 
In the course of our investigations, we discovered a subtle but severe mistake in the 


paper on the explicit determination of the outer (n—l)-radius of a regular simplex in 


E"^. Since this error completely invalidates the proof given there^, we give a description 
of that flaw, including some computer-algebraic calculations illustrating it. 

In that paper, the computation of the outer (n—l)-radius of a regular simplex is reduced 
to the analysis of the following optimization problem. 


(4.10) 


n+1 

min ^ uf 

i=l 

n+1 

s.t. Yl uf 
2=1 
n+1 

2 = 1 


1 , 

0 . 


2 


In a personal communication this has been confirmed by B. Weifibach. 
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For any local optimum (ui,..., Un+iY there exist Lagrange multipliers Ai, A 2 G M such 
that 

4'u^ -|- 2Ai'Uj -|- A 2 = 0, 1 ^ u -|~ 1' 

n+1 

(4.11) = 1, 

i=l 

n+1 

'Y^Ui = 0 . 

i=l 


Erroneously, in 0 


it is argued that symmetry arguments imply that A 2 = 0 in any 

shows 


solution. The following calculation in the computer algebra system Singular 
that for n = 3 this system has 26 solutions (counting multiplicity) over C. 


13 


ring R = 0, (ul,u2,u3,u4,lal,la2), (dp) 


ideal I = 

4*ul~3 + 2*lal*ul + la2, 

4*u2~3 + 2*lal*u2 + la2, 

4*u3~3 + 2*lal*u3 + la2, 

4*u4~3 + 2*lal*u4 + la2, 
ul~2 + u2''2 + u3''2 + u4~2 - 1, 
ul + u2 + u3 + u4; 


degree(std(I)); 

This program Erst dehnes a polynomial ring in the variables ui,... ,^ 4 , Ai, A 2 over a 
held of characteristic zero. We then use the degree command to compute the dimension 
and the degree of the ideal dehned by our equations. The output of that command is 


// codimension = 6 
// dimension = 0 
// degree = 26 


Hence, there are hnitely many solutions (since the dimension of the ideal is zero), and the 
degree of the ideal (the sum of the multiplicities of the solutions) is 26. 

18 of these solutions refer to the case A 2 = 0 (and those were the ones computed in 
Namely, if A 2 = 0 then the hrst row of (|4.11|) simplihes to 


Ui{2u‘^ + \i) = 0 , l<i<n+l. 


If we are only interested in the real solutions to this system, then setting Ai = —2A^ for 
some A > 0 gives 

Ui{uf — X^) = 0 , l<i<n+l. 
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Since the vector (mi, ..., Un+iY 
the solutions with A 2 = 0 are 


(0,..., 0)^ does not satisfy the second row in ([4.11D, 


Ui A, i G {^ 1 , • • • 5 ih\ 5 
Ui A, i G , ^2h} ) 

Mi = 0, f G {1,... ,n+1} \ {ii,. .. ,*2/1} 

for some h > 1, some set {zi,... ,Z 2 h} of pairwise different indices, and A = In 

the case n = 3, there are 12 possibilities to choose the indices and the signs for |h| = 1 and 
6 possibilities to choose the indices and the signs for \h\ = 2 , giving 18 solutions to ( |4.11|) . 

However, there are 8 additional solutions, which in fact are also real! Namely, these are 
the solutions 

1 7 1 

(ill,..., U 4 ) ^-^(1, 3,1,1) , Ai — , A 2 , 

1 7 1 

(mi, ..., U 4 ) 2 ^ 7 ^^ ^ 6 ' ^/3 ^ 


as well as the six distinct solutions obtained from these two by permuting the variables 
Ml,... ,M 4 . These additional solutions invalidate the subsequent arguments in 


31 


The omisssions get even worse in the higher-dimensional case. E.g., for n = 4, besides 
(D (D (D ( 2 ) = 20 -|- 30 = 50 solutions described in 
solutions: 


the 


31], we obtain the following 


(mi, ■ ■ ■ 1 ^ 5 ) 
(mi, ■ ■ ■ , U5 ) 
(mi, ■ ■ ■ , U5 ) 

(Mi, . . . ,M5)^ 


1 

Vso 

1 

1 

2^5 

1 

2V5 


(-2,-2,-2,3,3f, Ai = -^, X 2 = -^V^, 

15 75 

(2,2,2,-3,-3f, Ai = -^, = 

15 75 

1S 6 

( 15 — 4 , 1 , 1 , 1 )'^, 

13 fi 

(—15 4 , —1, —1, —1)^, ^ ~Yo ’ 


as well as those solutions obtained by permuting the variables. Altogether, we have 
10 -|- 10 -|- 5 -|- 5 = 30 solutions with A 2 7 ^ 0, and thus a total number of 80 solutions. 

Finally, we remark that the paper |^, which computes the outer (n—l)-radius of a 
regular simplex in odd dimension n, is correct (cf. also IQ). 
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